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Research status of density functional theory in corrosion of reactor alloy materials 
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Abstract After decades of research, the problems and behavior of the corrosion of reactor alloy materials are 


well known. However, some problems in corrosion of reactor materials have not been clarified, including the 


critical corrosion process of the reactor materials under operational conditions, the role of a single factor in 


the corrosion process, and the prediction of corrosion behaviors of new materials in extreme environments. 


The density functional theory, which is based on quantum mechanics, can accurately predict the motion 


process of atoms and the change in the relevant energy within a very short period. The density functional 


theory has become an important auxiliary method for investigating the corrosion process of reactor alloy 


materials in recent years and can help solve the above problems. In this paper, two parts are reviewed. The 
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first part introduces the density functional theory and includes the theoretical basis, development process, and 


mainstream computing software. In the second part, the research status of the density functional theory in the 


corrosion of reactor alloy materials is reviewed, including the adsorption, separation, combination, and 


internal diffusion of the reactor alloy material surfaces in the environments of water-cooled reactor, 


liquid-metal-cooled reactor, and molten salt reactor. 
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反应 堆 经 过 几 十 年 的 发 展 ， 已 经 衍生 出 多 种 推 型 ， 包 括 水 冷 堆 、 
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其 腐蚀 问题 主要 集中 在 表面 氧化 与 应 力 腐蚀 


开裂 方面 。 


目前 主要 采用 不 锈 钢 等 作为 包 壳 材料 ， 
第 一 壁 的 冷却 材料 ， 具 有 同样 的 问 
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高 温 氧 化 腐蚀 等 问题 。 
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可 以 精确 地 模拟 极 短 时 间 内 原子 的 


因素 所 起 的 作用 ; 


材料 以 及 极端 环境 下 材料 的 未 知 腐蚀 情况 ， 存 在 研究 难度 大 、 周 期 长 、 成 本 高 等 问题 。 


已 成 为 反应 推 合金 材料 腐蚀 过 程 和 
于 解决 上 述 问 题 . 尤 其 是 基于 量子 力学 的 密度 泛 函 理论 (Density Functional Theory, 
运动 过 程 与 能 量 的 变化 。 与 分 子 动力 学 、 相 场 、 有 限 元 分 析 等 其 


依然 存在 一 些 待 解决 的 问 
由 于 受到 观测 手段 的 限制 ， 实 验 上 难以 获得 腐蚀 中 最 本 质 的 反应 过 
二 ， 腐 蚀 过 程 往往 是 多 因素 耦合 的 复杂 行为 ， 难 


三 ， 对 于 新 


究 的 重要 辅 


他 模拟 方法 相 比 ，DFT 的 计算 结果 不 依赖 于 经 验 参数 ， 使 
金 材料 的 腐蚀 过 程 .实验 上 允 

这 为 研究 反应 堆 合金 材料 腐蚀 过 

究 通常 局 限 在 上 百 原子 的 空间 尺度 和 飞 秒 到 皮 秒 的 时 间 尺 度 内 。 


本 文 介绍 了 DFT 的 理论 基础 、 发 展 过 程 及 与 主流 计 犁 
反应 堆 合金 材料 腐蚀 模拟 中 的 应 用 情况 ， 包 括 水 冷 堆 、 液 态 金属 冷却 堆 以 及 炊 盐 地 


可 以 高 精度 的 模拟 不 同 条 件 下 反应 堆 合 


以 直接 观测 到 的 原子 尺度 腐蚀 现象 可 以 通过 DFT 模拟 直观 地 展示 出 来 ， 


程 及 机 理 提 供 了 重要 支持 手段 。 当 然 ， 作 为 高 精度 的 代价 ，DFT f 
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并 总 结 了 DET Wiz 


在 多 


密度 泛 函 理论 (DFT) 是 为 了 求解 多 粒子 薛 定 齐 方 程 而 提出 的 ， 如 下 所 示 : 


Hr, R) = E"V(r , R) 
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能 与 电子 -原子 核 相互 作 


解 ，Born 和 


EE 子 的 坐标 ;，R 代表 所 有 原子 核 的 坐标 。 方 程 包括 
j。 众 所 周知 ， 电 子 与 原子 核 的 质量 与 速度 悬殊 ， 为 了 简化 薛 定 齐 方 程 的 求 
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外 子 的 动能 、 原 子 核 的 动 


Oppenheimer 提出 了 绝热 近似 帆 。 通 常情 况 下， 电子 在 高 速 地 运动 ， 原 子 核 的 运动 远 远 
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跟 不 上 电子 的 运动 。 因 此 ， 当 我 们 考虑 电子 的 运动 时 ， 可 以 认为 原子 核 产 生 的 势 场 是 无 限 组 变 的 ， 
从 而 可 以 分 别 计算 电子 和 原子 核 的 运动 。 

1964 年 ，Hohenberg 和 KohnP 3E tH] Hohenberg-Kohn 定理 是 DFT 的 基础 。Hohenberg-Kohn 定 
理 的 基础 概念 是 Thomas 和 Fermi 提出 的 粒子 密度 函数 BI， 它 可 以 分 成 两 个 基本 定理 : 一 个 是 多 粒 
子 薛 定 齐 方 程 中 的 基态 能 量 ， 由 粒子 密度 唯一 决定 ;， 另 一 个 是 粒子 数 恒定 下 粒子 密度 使 总 泛 函 的 能 
量 最 小 ， 也 就 是 基态 能 量 。Hohenberg-Kohn 定理 中 基态 能 量 泛 函 可 以 表达 为 : 

E[(Z] = Exow [UP] + Exc [UP] Q) 

SUH: Exc[ 中 是 交换 关联 能 泛 函 : Erow [四 分 为 以 下 4 部分: 


h2 
Exsown[{Pi}]=- —— YW; V?W , d ?r +fV (r )n (r )dr + 


e 26n (r )n (r ^) 43 3 & 
j ff TS d?r d?r + Eion (3) 


从 左 到 右 分 别 为 电子 的 动能 、 电 子 与 原子 核 的 相互 作用 、 电 子 对 的 相互 作用 、 原 子 核对 的 相互 作用 。 
1965 ££, Kohn 和 Sham 等 人 外 为 了 求解 基态 能 量 泛 函 ， 提 出 了 Kohn-Sham 方程 ， 如 下 所 示 : 


pr)- SAI, (r dP (4) 
CV Verl P Eh VU. =E i V, (r ) (5) 
r ` 6 X 
Valp(r)- v (r )+ fdr v — (6) 


AO 右 侧 分 别 是 电子 动能 、 库 伦 势 和 交换 关联 势 。 上 面 的 (4)、 C50 I (6) 三 式 一 起 称 为 Kohn-Sham 
方程 。 求 解 时 ， 设 定 一 个 初始 电荷 密度 后 ， 通 过 方程 《4) SR, (7 ”)， 然 后 带 


入 方程 (5) 计算 电荷 密度 ， 并 与 初始 电荷 密度 进行 比较 。 如 果 两 个 密度 相同 ， 这 就 是 基态 电子 密度 ， 
可 以 用 来 计算 总 能 量 。 但 在 实际 中 ， 相 互 作用 电子 的 交换 相关 能 泛 函 Exc[ p] 5 无 法 精确 得 到 。 为 了 
求解 Kohn-Sham 方程 ， 必 须 找到 交换 关联 能 泛 函 。Kohn 和 Sham 等 人 四 提出 了 均匀 电子 气 的 概念 ， 
并 通过 局 部 密度 近似 (LDA: Local Density Approximation) PJZ REKREI Ko LDA 仅 适 用 于 均匀 
电子 体系 ， 为 了 应 用 于 非 均 匀 电 子 体系 ， 在 LDA 基础 上 加 入 电子 密度 梯度 ， 得 到 了 广义 梯度 近似 
(Generalized Gradient Approximation, GGA)。 其 表达 式 很 多 ， 例 如 BLYP、BP9%6、PW91 和 PBEUIAE, 
由 于 LDA/GGA 泛 函 高 估 了 电子 离 域 ， 难 以 描述 过 渡 金 属 如 铬 等 含有 d. f 轨道 电 子 的 强 关 联 体系 ， 
因此 将 电子 分 为 两 个 子 系统 : 对 于 s EXE p 轨道 电子 ,仍然 用 LDA/GGA 泛 函 来 描述 ; 对 于 局 域 d 或 f 
轨道 的 电子 ， 通 过 引入 Hubbard 项 ， 增 加 在 位 库伦 修正 (on-site coulomb correction) 来 处 理 。 此 方 
法 称 为 DFT+U 方法 外 。 在 固体 物理 学 计算 中 ， 价 电子 决定 原子 的 电子 结构 性 质 。 因 此 ， 为 了 节约 计 
算 资 源 ， 采 用 族 势 方法 对 单 电子 方程 进行 改造 。 即 价 电子 本 征 函 数 如 实 近 似 ， 非 价 电子 本 征 函 数 平 
滑 处 理 。 常 用 的 帮 势 为 : STERA (Norm Conserving Pseudopotential, NCPP )、 超 软 帮 势 (Ultrasoft 
Pseudopotential, USPP〉BP71、 投 影 级 加 平面 波 〈Projector-augmented Wave, PAW) 0 等 。 

在 使 用 DET 研究 反应 堆 合 金 材料 的 腐蚀 问题 中 , 主要 关注 的 是 在 不 同 条 件 下 反应 堆 合 金 材料 腐 
蚀 行 为 中 的 体系 能 量变 化 、 化 学 键 重 构 和 电荷 转移 等 ; 研究 对 象 包括 反应 堆 合 金 材 料 的 基体 、 唱 面 、 
晶 界 和 界面 等 ， 研 究 的 化 学 反应 包括 反应 介质 在 表面 的 吸附 及 在 基体 与 晶 界 的 扩散 等 ， 具 体外 界 条 
件 因素 包括 应 力 、 温 度 、 辐 照 等 等 。 其 中 ， 从 头 算 分 子 动力 学 (Ab Initio Molecular Dynamics, AIMD ) 
可 以 实现 有 限 温度 下 腐蚀 行为 的 研究 。 在 辐 照 对 腐蚀 影响 这 方面 ， 由 于 DFT 研究 尺度 过 于 微观 ， 该 
领域 的 DFT 研究 较为 有 限 。 目 前 主要 研究 各 种 点 缺陷 对 腐蚀 的 影响 ， 其 他 辐 照 现 象 〈 例 如 : 辐 照 肿 
胀 、 辐 照 硬 化 、 辐 照 蠕 变 和 辐 照 脆 化 等 ) 对 腐蚀 行为 影响 的 DFT 计算 未 见报 道 ， 这 也 是 今后 使 用 理 
论 模拟 方法 研究 反应 堆 合 金 材料 腐蚀 的 发 展 方向 。 除 了 明确 看 到 的 原子 位 移 之 外 ，DFT 计算 可 以 得 
到 的 结果 包括 表面 能 、 吸 附 能 、 功 函 等 能 量 信息 和 态 密度 (Density of States, DOS )、 能 带 、 差 分 电 
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另外 ， 通 过 AMD 与 过 渡 态 搜索 还 可 以 寻找 反应 中 间 态 结构 3 


随 着 计算 机 技术 的 发 展 和 计生 


方法 的 改进 ，DFT 应 月 


ts 


相关 原子 


有 取得 了 长 足 进步 ， 相 关 文 章 发 表 数 快速 增 
Ko 与 此 同时 ，VASP (Vienna Ab initio Simulation Package). Materials Studio. Gaussian 等 商业 软 伯 


FE 和 


Quantum Espresso 、SIESTA(Spanish Initiative for Electronic Simulations with thousand of Atoms) 、 


ABINIT 等 开源 软 伯 
ATE, RS 
间 5 #4 

Studio fll SIESTA. 
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F 也 被 研究 人 员 用 于 材料 性 能 的 下 
LE 论 基 础 主要 分 为 平面 波 、 原 子 轨 道 线性 


合 、 线 性 缀 加 3 


三 种 经 常 使 月 


究 模 拟 。 表 102 给 出 了 目前 基于 DFT 的 主流 计 
FHE. B 
HA) DFT 计算 软件 : 


1 展示 了 近 20 年 
VASP. Materials 


VASPUS HUE SC HS TETUER AA LDA 和 GGA 两 种 交换 关联 泛 函 来 描述 电子 - 核 相 互 作 用 的 计算 
软件 。VASP 可 以 计算 固体 表面 和 液体 系统 中 原子 、 分 子 、 团 徐 、 纳 米 结构 的 结构 参数 、 键 参数 、 力 
学 性 能 、 光 学 和 磁性 性 能 ， 在 反应 堆 合 金 材料 腐蚀 计算 中 应 用 最 广泛 。Materials Studio 是 兼 具 建 横 
和 计算 的 综合 性 平台 ，CASTEP 和 DMol3 是 基于 DFT 模拟 计算 的 模块 (9。CASTEP 采用 平面 波 大 
势 和 LDA M GGA 两 种 交换 关联 泛 函 进行 计算 ， 可 用 于 研究 固体 材料 的 能 带 结 构 、 态 密度 、 电 和 荷 分 
布 、 表 面 化 学 等 。DMol3 基于 量子 化 学 ， 主 要 用 于 研究 表面 和 固体 材料 的 均 相 催化 、 多 相 催 化 、 分 
FRM. SIESTA 是 一 个 基于 DFT 的 开源 软件 包 061, 它 采 用 标准 的 Kohn-Sham 方法 , 使 用 LDA/GGA 
来 计算 体系 能 量 与 Mulliken 电 蓓 分 布 、 原 子 轨道 、 电 蓓 密度 、 能 带 和 态 密度 等 。 

表 1 以 DFT 为 基础 的 软件 简介 
Table 1 Introduction to software based on DFT 
软件 名 称 使 用 许可 
Name License 
平面 波 为 基础 
Plane-waves basis sets 
TES H-a 
VASP 商业 化 软 人 
Commercial software ? 
Quantum Espresso 基于 GPL 协议 使 用 
p Based on GPL 
TET H- b 
CASTEP 商业 化 软 人 
Commercial software? 
ABNT 基于 GPL 协议 使 用 
Based on GPL 
POR 基于 GPL 协议 使 用 
Based on GPL 
Be 
CPMD 免费 使 用 
Free 
ONETEP 商业 化 软件 
Commercial software 
BigDET 基于 GPL 协议 使 用 
Based on GPL 
原子 轨道 线性 组 合 为 基础 
Atom-centered basis sets 
商业 化 软件 
Gaussian 
Commercial software 
GAMESS 免费 使 用 
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Free 
商业 化 软件 
Commercial software 
免费 使 用 。 
Free* 
商业 化 软件 
Commercial software 
免费 使 用 
Free 
商业 化 软件 5 
Commercial software 
商业 化 软件 
Commercial software 
商业 化 软件 


Commercial software 


Molpro 


SIESTA 


Turbomole 


ORCA 


CRYSTAL 


Q-Chem 


FHI-aims 


实 空间 格 点 为 基础 
Real-space grids 


基于 GPL 协议 使 用 
Octopus 
Based on GPL 
Gd 基于 GPL 协议 使 用 
Based on GPL 
线性 级 加 平面 波 : 
WIEN2k 商业 化 软件 
Commercial software 
基于 GPL 协议 使 用 
exciting 
Based on GPL 
Be 
EEBUR 免费 使 用 
Free 


TE: a- 澳 大 利 亚 学 术 机 构 免 费 ，b- 英 国学 术 机 构 免 费 ，c- 学 术 机 构 免费 ，d-CP2K 混合 了 平面 波 与 原 
子 轨道 线性 组 合 方法 ，e-GPAW 也 可 以 使 用 平面 波 或 原子 轨道 线性 组 合 方法 


Notes: a-Free for academic institutions in Austria; b- Free for academic institutions in UK; c- For 


academics; d- CP2K employs mixed plane-waves and atom-centered basis sets; e-GPAW can also employ 


plane-waves or atom-centered basis sets 
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图 1 主流 DFT 计算 软件 的 文章 发 表 数 〈 数 据 来 自 Google Scholar) 
Fig.1 Number of articles published by mainstream DFT computing software (The data are obtained from Google Scholar) 


反应 堆 合金 材料 在 水 冷 堆 环境 下 腐蚀 的 DFT 模拟 计算 


水 冷 堆 是 目前 在 和 运行 反应 堆 中 占 比 最 大 的 一 种 堆 型 。 对 水 冷 动 力 堆 主 要 结构 材料 锋 合 金 、 不 锈 
钢 及 钊 基 合 金 在 水 环境 下 的 腐蚀 行为 实验 研究 时 间 较 长 〈 已 超过 50 年 )， 腐 蚀 行为 较为 明确 ， 相 关 
腐蚀 试验 研究 文献 数 以 千 计 ， 但 是 对 其 腐蚀 机 理 的 理解 依然 不 够 明确 。 近 些 年 已 有 学 者 开展 水 冷 堆 
环境 下 反应 堆 合 金 材料 腐蚀 机 理 的 DET 研究 ， 然 而 ， 相 比 实验 工作 而 言 ， 相 关 DET 计算 的 研究 工 
作 仍 然 非常 有 限 。 表 2 汇总 了 水 冷 堆 环 境 下 反应 堆 合 金 材料 腐蚀 DFT 计算 的 文献 情况 。 其 中 ， 大 多 
数 的 腐蚀 计算 工作 围绕 钳 的 腐蚀 问题 开展 ， 另 有 少量 FeCrA1、 不 锈 钢 的 腐蚀 计算 。 
表 2 水 冷 堆 环境 下 反应 堆 合金 材料 腐蚀 DFT 计算 文献 汇总 
Table 2 Literature on DFT calculations of corrosion of reactor alloy material in water-cooled reactor environment 
研究 机 构 第 一 作者 主要 内 容 年 份 
Research institution First author Research content Year 
Hunan university Z. Liu 氧 在 Zr(1011) 表 面 的 吸附 行为 0 Adsorption 2021 
behavior of oxygen on Zr (101 1) surface 
Capital Normal University F. Wang 氧 在 Zr(0001) duré ARPES) Adsorption 2008 
behavior of oxygen on Zr (0001) surface 
Shanghai University H. Zhang UTE RAR AT PE v TAL ES] V ERR UI 
Adsorption of oxygen on zirconium surface 2018 
with low Miller index 
University of Oxford R.J. Nicholls — $555/8/ SU FR LE Ab R E UP 
的 结构 PO Structure 2013 
of the zirconium oxide phase found at the 
zirconium metal/oxide interface 
Kim Il Sung University J.Y.Jong BERTA ren ae P Pe AE AK SE AY 
吸附 作用 
Adsorption of zirconium to hydrogen 2023 
produced during corrosion and hydrogen 
co-existing in water 
Idaho National Laboratory M.V. Glazoff FEE um ic io AY EA 
The role of zirconia in hydrogen absorption 2014 
and diffusion 
Université Paris-Saclay E. Haurat PUR SUL AHI BUT AP] Hydrogen 2022 
diffusion behavior in monoclinic zirconia 
Hunan university M. Feng ZNJI SE THEXETA SS (hcp-Zr) rP Z8 -AHE 
{EAP Impurity- 2020 
Hydrogen interaction in hcp-Zr 
Xi'an Jiaotong University S.M.Liu ATE 8 XE PR Pie PESE REUS! The 
process of precipitation of hydride in 2022 
zirconium matrix 
Institute of Applied Physics and X.Y. Zhu ZrHx 的 热力 学 性 质 (x = 0.5,1,1.5,2)06 2018 


Computational Mathematics 
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Thermodynamic properties of ZrHx 


(x 


University of Tokyo 


Malmo University 


Universite' de Lille 


Imperial College London 


State Power Investment 
Corporation Research Institute 


Los Alamos National Lab 


Los Alamos National Lab 


North Carolina State University 


City University of Hong Kong 


Sun Yat-Sen University 


Sun Yat-Sen University 
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Y. Udagawa 


P.A.Olsson 


A. Legris 


V. Podgurschi 


R. Tu 


M. L. Rossi 


M. L. Rossi 


Rák 


J. Wei 


XJ. Li 


XJ. Li 


= 0.5,1,1.5,2) 


45-55 t Le Br 22 VE BE AT EELT Comparison 
of fracture properties of zirconium and 
hydride 

SURGE T DN J ERER SLS Effect 
of hydrogen filling vacancy on mechanical 


properties of zirconium 


碘 在 错 基 面 与 棱柱 面 的 吸附 行为 29 


Adsorption behavior of iodine on zirconium 


base and prismatic surface 


施加 静水 压力 (-2% ~ +3%) 条 件 下 ， 碘 和 氧 
与 错 基 体 和 表面 的 相互 作用 B09] Interaction 


of iodine and oxygen with zirconium substrate 


and surface under hydrostatic pressure (—296 ~ 
+3%) 


EE Ee RUE IT ME. A RAR 
频率 BU Formation 


energy, diffusion barrier and vibration 


frequency of complex iodine defects in 


Zirconium 
不 同 的 PBE EAER E P E r 
误差 的 对 比 B2 Error 


comparison of different PBE 
pseudopotentivities in calculation of various 


zirconium iodide products 


ISCC HEI W- SC- 75 AY Be I RECS 
Iodine-oxygen-zirconium reaction process 
during ISCC 


奥 氏 体 不 锈 钢 (304\316) 的 表面 能 B4 Surface 


energy of austenitic stainless steels (304\316) 


水 在 FeCrA1(110) 表 面 吸附 的 的 结构 B3] 
Atomic structure of water adsorption process 
on FeCrAI (110) surface 


AE Fe(100) 和 FeCrAI (100) 表 面 的 吸附 B9 
Hydrogen adsorption on Fe (100) and FeCrAI 
(100) surface 


Fe, Cr. Al 和 FeCrAl 晶体 (110) 表 面 上 氧 原 
子 的 吸附 行为 B7 Adsorption 
behavior of hydrogen atoms on (110) surfaces 
of Fe, Cr, Al and FeCrAI 


2010 


2015 


2005 


2022 


2018 


2013 


2015 


2017 


2019 


2021 


2022 


2021 
Anhui University Z. Liu Fe £ 3 «110» (111) 和 5 <001> (310) HFEF 


氧 的 特性 B39 Properties 
of oxygen in grain boundaries of Fe X3 <110> 
(111) and X5 «001» (310) 


13 ” 钳 合 金 表面 氧化 问题 的 DFT 计算 


独 合金 在 水 冷 堆 环境 的 表面 氧化 过 程 为 钳 与 水 反应 生成 氧化 错 ， 同 时 解 离 出 氧 。DEFT 计算 内 容 
包括 氧 原 子 和 氧 分 子 在 钳 表 面 的 吸附 、 氧 在 钳 表 面 的 扩散 、 猪 氧化 物 薄膜 的 晶体 形 貌 和 生长 速率 、 
告 氧化 物 的 相 变 等 。 

在 错 合 金 腐蚀 过 程 中, 氧 原子 首先 吸附 于 钳 合 金 表面 .在 钳 合 金 (0001) 表 面 吸附 位 点 如 图 2 所 示 ， 
绿色 原子 为 表面 错 原 子 ， 蓝 色 原 子 为 次 表面 猪 原子。 红色、 黄色 、 粉 色 、 青 色 原 子 为 吸附 在 表面 不 
同位 点 的 氧 原子 。 红 色 原 子 位 于 表层 钳 原 子 正 上 方 ， 该 吸附 位 置 称 为 TOP 位 ; 粉色 原子 位 于 表层 三 
个 独 原 子 中 心 ， 该 吸附 位 置 称 为 FCC 位 ; 黄色 原子 位 于 表层 两 个 铬 原子 中 心 ， 该 吸附 位 置 称 为 BRI 
位 ; 青色 原子 位 于 表层 三 个 独 原 子 中心 与 次 表层 钳 原 子 正 上 方 ， 此 吸附 位 置 为 HCP 结构 特有 的 吸附 
位 点 ， 因 此 称 为 HCP 位 。Liu07 等 人 使 用 了 GGA-PBE 泛 函 进行 DFT 模拟 ,研究 了 氧 在 铬 (1011) 表 
面 的 吸附 与 扩散 行为 。 结 果 表 明 ， 钳 的 4d 价 电子 倾向 于 迁移 到 被 吸附 0; 分 子 的 反 键 轨道 上 ， 导 致 
O-O 键 断 裂 并 释放 能 量 。Wangl'5 等 人 使 用 GGA-PW91 泛 函 进行 DFT 模拟 在 计算 氧 原子 在 钳 (0001) 
唱 面 的 吸附 过 程 中 , 发 现 氧 覆盖 率 为 0.2SML 情况 下 FCC 位 点 是 最 稳定 位 点 , 随 着 氧 覆 盖 率 从 0.2SMEL 
增加 到 1.0ML, 每 个 氧 原子 倾向 于 从 上 往 下 依次 占据 各 层 的 FCC 位 点 。Zhang09 等 人 使 用 GGA-PBE 
i BERT DFT 模拟 研究 了 氧 在 三 个 低 Miller 指数 错 表 面 的 吸附 规律 。 结 果 表 明 ， 在 氧 迅 速 减少 的 条 
件 下 和 氧 优先 吸附 在 钳 基 体 表 面 (0001)， 而 在 氧 缓慢 减少 和 氧 富 集 的 条 件 下 氧 优 先 吸 附 于 在 棱柱 表面 
(1120) 和 (1010)。 

氧 在 铬 表面 吸附 扩散 后 ,与 钳 快 速 结合 形成 致密 ZrO» 氧化 层 导致 氧化 速率 降低 ， 再 之 后 氧化 速 
率 增加 。NichollsPR9I 等 人 认为 这 可 能 与 错 亚 氧化 物 的 形成 有 有关， 结合 DFT 和 高 分 辨 率 电 子 显微镜 ， 
在 钳 金 属 /氧化 物 界面 处 发 现 了 ZrO 亚 氧 化 物 的 结构 亚 氧化 物 计算 得 到 的 峰 位 置 与 电子 能 量 损 失 谱 
测 出 的 峰 位 置 接近 。 亚 氧化 物 能 分 解 为 较为 稳定 的 Zro 和 ZrO。 随 着 腐蚀 的 进一步 发 展 ， 金 属 / 氧 
化 物 界面 〈 发 生 电 化 学 反应 的 点 ) 逐渐 向 金属 基体 方向 移动 ， 远 离 新 形成 的 氧化 物 。 随 着 远离 该 界 
面 的 距离 增加 ， 氧 化 层 中 的 压 应 力 显 着 减 小 ， 其 中 应 力 不 再 足以 稳定 氧化 层 结构 后 ， 导 致 氧化 层 开 
裂 ， 从 而 氧化 速率 迅速 增加 。 


图 2 钳 (0001) 表 面 常见 吸附 位 点 (彩色 见 网 络 版 
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Fig.2 Common adsorption sites on Zr (0001) surface (color online) 
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Zr 表面 (0001) 的 吸附 现象 。 发 现 两 者 最 有 利 的 吸附 位 点 都 是 hop 位 点 。 当 
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1.5” 钳 合金 ISCC 问题 的 DFT 计算 
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会 与 错 反 应 生成 Zril, Zrl3 和 Zr[4， 


能 和 不 稳定 堆 
能 小 23%， 容 易 产 生 断 裂 面 ， 


以 进行 位 错 运动 ， 因 此 氢化 锋 容 易 脆 化 。OlssonPs 等 人 使 


。 由于】 


包 壳 的 性 能 。 


IZ 


然 料 裂变 过 程 


E 高 燃 耗 或 热 应 力 导 致 芯 块 开裂 的 情况 下 ， 该 ; 


能 的 影响 ， 发 现 空位 中 的 氧 降 
层 错 不 稳定 ， 从 而 降低 了 塑性 。 空 位 的 引入 降低 了 断裂 的 峰值 应 力 ， 而 空位 被 


Ex 


EST 


中 Zr 在 200-500 %C 下 具有 挥发 性 , 在 1100-1 
的 地 方 并 通过 化 学 平衡 释放 分 子玉 


碘 不 断 
过 程 会 加 


人 


Legris’? 


| GGA Z EAT f DFT 模拟 , WEF T E 


与 外 的 相互 作 上 月 


SRI RE | 


5n 


已 


A Ta EA 


用 GGA-PBE 32 AE 


HT WER (080% 以 上 )， 基 二 


WEH 


H 


“HABA 3.2 10 5 em? s! . Podgurschi 


20 等 人 认为 氧 也 影响 


ÍT DFT RM, META 


TS ESTER ELTE H 
AA 


H 


成 能 对 应 变相 对 不 敏感 。 随 着 


结果 表明 ,fn 


p 


E 


S 


f SUM 


间 质 的 保护 作 
泛 函 在 300~2 000 K 的 温 
扩散 势 垒 和 振动 频率 ， 有 具体 缺陷 如 图 


Jo [x] 时 4A 
度 范 围 内 , 采 月 
4 所 示 


E Jj JI 
i ds Per /\ TETA TEI BR BE EZ aR Re BE DL RE E BA CC 


RTRA SERRER. Tuen 
有 AIMD FEIRA KHH T BAH OO 


, FIN 


Hl 


日 ,结果 表明 看 
能 降低 比 棱柱 面 要 高 ， 这 与 实验 结果 一 致 。 
J ISCC 的 发 展 , 并 使 


AI 人 YE 


吸附 导致 


H 


[流体 静 压 (-2% 至 +3%) 条 


F 下 裂纹 


mà 


IT 


端 处 碘 - 氧 - 


UU. T 
的 增加 ， 碘 间 质 的 形成 能 从 1.03 eV 增加 到 8. 


jt bud és 
61 eV， 证 明 


H 


di 
pu 


1 GGA-PBE 


HR 


E £h 


原子 为 碘 原 子 取 代 了 铬 原子 并 将 


C4 


置 ， 


你 为 置换 填 


隙 型 缺陷 粉色 原子 为 


fi 
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原子 占据 了 间隙 位 置 ， 称 为 
构成 了 研究 辐 照 点 缺陷 对 反应 夫 
间隙 扩散 和 带 有 空位 的 置换 型 


间隙 型 缺陷 ; 


取代 了 钳 原子 的 位 置 ， 称 为 
线 代 表 锋 原子 确实 ， 称 为 空位 型 缺 


I: 


会 


ELX 


P e 


自 的 基本 构 型 。 发 现 300~2 000 K 的 


1 


[名 


陷 的 形成 能 、 


Fi 


移动 至 间 隐 位 


(44 
EE 


黄色 原子 为 
。 四 种 缺陷 
围 内 ， 碘 的 


扩散 机 制 ， 间 隙 扩散 最 快 ， 空 位 会 减 绥 


i 


原子 的 扩散 。 


图 4 钳 晶 体 中 常见 缺陷 。 绿 色 原 子 为 表面 错 原 子 ， 蓝 色 原 子 为 次 


表面 钳 原 子 ， 虚 线 为 表面 钳 原 子 原 本 在 的 位 置 。 紫 


fi, RE. EIR TAA FA oR AE 


原子 (彩色 见 网 络 版 ) 


Fig. 4 Common adsorption sites on Zr (0001) surface. Green and blue atoms are the surface and sub-surface zirconium atoms, respectively. Purple, yellow, pink 


atoms are adsorbed iodine atoms that cause different defects 

研究 ISCC PLAY, Rossi B3 等 人 考虑 了 Zr OPE T 

的 PBE 硒 势 、 范 德 华 力 色散 校正 的 PBE RABE HTL, ACE EA 
到 的 体积 模 量 与 真实 情况 最 为 接近 。 之 后 ，RossiB3] 等 人 推断 蝇 


(color online) 
E, 将 未 矫正 的 PBE 厦 势 、 半 经 验 色散 校正 
范 德 华 力 色 散 校正 的 PBE REST 
H ISCC 期 间 可 能 的 反应 过 程 ， 由 于 钳 


表面 存在 保护 性 氧化 物 ， 碘 无 法 直接 蔡 换 表面 的 氧 。 然 而 ， 存 在 缺陷 的 情况 下 碘 可 以 沿 外 基体 表面 


形成 ZrIx 相 的 混合 物 ， 形 成 氧化 物 - 碘 化 物 系 统 。 随 着 环 向 应 力 的 增加 ， 富 含 


生 点 蚀 ， 最 终 成 为 应 力 腐蚀 开裂 的 起 始点 。 


L6 其 他 反应 堆 合金 材料 水 腐蚀 的 DFT 计算 


RiikB4 等 人 使 用 GGA-PBE 泛 函 进 行 DFT 模拟 ,研究 了 两 
(100)、(110) 表 面 能 ， 发 现 奥 氏 体 不 锈 钢 存 在 各 向 异性 ， 排 列 最 紧密 且 最 稳定 的 是 (111)， 其 次 是 (100) 


和 (110) 表 面 。 作 者 建立 了 基于 温度 、pH 值 与 水 溶液 中 离子 (Fe 六 、 
表面 稳定 性 的 模型 。 发 现 (111) 面 的 表面 稳定 性 随 温 度 和 pH 值 的 升 高 
少量 增加 。WeiB3 等 人 使 用 GGA-PBE 泛 函 进行 DFT 模拟 , 研究 了 水 在 FeCrAI (110) 表 面 吸附 的 各 


Fet, NIH Cr? ik AIT AS A 
[降低 ; 随 水 溶液 离子 浓度 增加 


构 型 ， 结 果 表 明 水 在 FeCrAl (110) 表 面 的 Top 位 吸附 能 最 


等 人 使 用 GGA 研究 了 氧 原 子 和 分 子 在 铁 (100) 和 FeCrAI(100) 


大 ， 元 素 


的 唱 界 容易 断裂 并 发 


! 奥 氏 体 不 锈 钢 (304 和 316 型 ) 的 (111)、 


1 铬 元 素 与 水 的 结合 能 更 大 。LiBg 
表面 的 吸附 ， 发 现 铝 原 子 本 身 与 氧 原子 


没有 相互 作用 ， 铬 原子 旁边 的 铝 原子 会 增强 铬 原子 与 氨 的 相互 作用 。 当 铬 原子 在 铁 原 子 劳 边 时 ， 这 


两 个 原子 ( 铁 - 铬 桥 位 ) 可 以 减少 与 氨 的 相互 作用 。LiB? 等 人 后 续 使 用 GGA-PBE 泛 函 进行 DFT 模拟 ， 


研究 了 铁 、 铬 、 铝 和 FeCrAl 四 种 体 心 立方 晶体 (110) 表 面 上 氧 原子 
氢 的 吸附 有 较 强 的 影响 ， 铝 原子 主要 通过 影响 铬 与 氢 的 相互 作用 来 影响 吸附 。Liu 时 等 人 使 
GGA-PBE 泛 函 进 行 DET 模拟 , 研究 了 体 心 立方 铁 晶 界 23<110>(111) 和 5<001>(310) 中 间 阶 氧 的 溶 


和 扩散 特性 ， 间 隙 氧 原子 更 倾向 于 溶解 在 体 心 立 方 铁 中 ， 
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Id 


E dn 


E58 


Hi 
中 


的 吸附 行为 。 结 果 表 明 ， 铬 和 铝 对 


BS i 


局 析 从 而 氧化 唱 界 。 


反应 堆 合金 材料 在 液态 金属 冷却 堆 环境 下 腐蚀 的 DFT 模拟 计算 


液态 金属 具有 导热 系数 高 、 低 熔点 、 高 沸点 、 低 蒸汽 压 等 优点 ， 也 是 反应 堆 


如 ， 液 态 钠 
用 作 铅 基 快 堆 的 冷却 剂 09， 


的 理想 冷却 剂 。 例 
j 作 钠 冷 快 堆 的 冷却 剂 B9， 液 态 铅 或 液态 铅 锐 共 晶 合 金 〈Lead-Bismuth Eutectic, LBE) 


WAS BEA (EAR AR HEA A AN FS ISS. 液态 金属 快 中 子 反 应 堆 中 ， 


燃料 包 壳 及 结构 材料 多 选用 不 锈 钢 等 铁 基 合金 。 在 材料 相 容 性 方面 ， 液 态 金 属 之 间 存 在 着 不 同 的 特 


性 名。 例如 ， 在 液态 馈 和 LBE 中 不 锈 钢 很 容易 发 生 溶解 腐蚀 ， 尤 其 是 其 
液态 钠 《〈 氧 含量 控制 在 10 ugg! 以下) 与 钢 的 相 容 性 较 好 比 


金 材料 腐蚀 的 部 分 DFT 文献 。 


。 表 3 汇总 


的 钊 元 素 中 ]， 但 是 高 纯度 


了 液态 金属 环境 下 反应 堆 合 


表 3 液态 金属 环境 下 反应 堆 合 金 材料 腐蚀 DFT 计算 文献 汇总 
Table 3 Literature on DFT calculations of corrosion of reactor alloy material in liquid-metal environment 
研究 机 构 第 一 作者 主要 内 容 年 份 
Research institution First author Research content Year 
Seoul National University J. Gil 氧 在 液态 Na "PES TR ARM! 2020 
Enthalpy of dissolution of oxygen in liquid 
sodium 
Japan Atomic Energy Agency S. Kikuchi ”水 分 子 在 液态 Na 表面 的 反应 过 程 [0 2012 
Reaction of water molecules on the Surface of 
liquid sodium 
POSTECH, Republic of Korea S. J. Kim 钠 在 Ti(0001) 表 面 的 吸附 的 Adsorption 2015 
of sodium on Ti (0001) surface 
Ohio State University X.Li 液体 钠 中 挫 杂 铀 系 元 素 的 影响 [9 2017 
Effect of lanthanide doped in liquid sodium 
Ohio State University A. Samin 1000 K 温度 下 液态 钠 中 锅 的 性 质 的 2015 
The properties of cerium in liquid sodium at 
1000 K 
Sichuan University X. Long 铅 和 锐 原 子 在 Fe(111) 表 面 的 相互 作用 60 2022 
Interaction of Lead and bismuth on Fe (111) 
surface 
Institute of Solid State C.Song LBE 与 Fe(100) 表 面 的 相互 作用 5 2014 
Physies, Chinese Academy of Interaction of LBE on Fe(100) surface 
Sciences 
Sichuan University Y.FLi LBE 原子 与 Fe(11 D/EDBSZE IAA EEE?) 2023 
Interaction of LBE atoms with rough Fe(111) 
surface 
Sun Yat-Sen University T. Liu LBE 在 铁 素 体 钢 晶 界 处 的 腐蚀 机 理 及 合金 元 2022 
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素 的 影响 3 

Corrosion mechanism of LBE at the grain 
boundary of ferritic steel and influence of 
alloying elements 


Sun Yat-Sen University T. Liu 
Sichuan University R.Zhou 
Chinese Academy of Sciences W. Ding 
Chinese Academy of Sciences Y. Xu 
Tiangong University Y.Li 


1.7 钠 冷 环境 中 钠 的 DFT 计算 
关于 钠 冷 快 堆 钠 与 不 锈 钢 相 容 性 问题 ， 


冷 快 堆 的 工程 化 设计 和 建造 ， 开 展 了 大 量 的 验证 


Fe 和 Fe304 表面 对 LBE 腐蚀 的 敏感 性 [5 


Corrosion sensitivity of Fe and Fe 
to LBE 


合 氧 LBE 与 含 合金 元 素 的 Fe 表面 的 腐蚀 行为 


[55] 


behavior of oxygen-containing LBE and Fe 


surface with alloying element 


合金 元 素 与 氧 在 钢 -LBE 界面 上 的 行为 59 


Behavior of alloying elements and oxygen at the 


steel-LBE interface 


锂 在 铁 表面 的 吸附 和 铁 向 液态 外 


Adsorption of lithium on the surfa 


the escape of iron to liquid lithium 


锂 在 bcc-Fe 中 扩散 的 性 质 68] 


Properties of lithium diffusion in bcc-Fe 


国内 外 学 者 展开 了 较为 


题 鲜 见 文献 报道 。 
反应 的 模拟 研究 。 


— 


且 是 ， 针 对 钠 冷 快 堆 的 运行 与 安全 ， 确 


泛 的 试验 在 
性 考验 ， 但 是 利用 DFT 方法 预测 名 
有 学 者 利用 DFT FET $ 


2023 
304 surfaces 
2023 
Corrosion 
2019 
逃逸 过 程 6] 2019 
ce of iron and 
2009 
究 ， 尤 其 是 针对 钠 
中 材料 的 腐蚀 问 
的 性 质 以 及 钠 水 


Gil SYEH AIMD 对 钠 物 性 进行 了 模拟 研究 , 得 到 在 600 K 和 1000K 温度 下 氧 在 液态 钠 中 的 溶 


MK, 


kJ-mol"! 和 38.2 kJ-mol! 等 计算 结果 ,通过 与 实验 数据 比 对 , 证 明 计 算 值 与 实验 值 


从 而 验证 


了 该 方法 的 有 效 性 。 


FERNY p 
发 现 反 应 过 程 中 水 分 子 所 具有 的 平均 动 
结果 表明 : 水 在 液体 钠 表面 
子 朝向 液体 钠 表 面 
SWR， 并 计算 了 钠 在 Ti(0001) 面 的 吸附 能 ， 
面 易 相互 结合 ， 导致 在 钛 纳米 粒子 表面 上 吸 


钠 一 旦 发 生 泄露 , 与 空气 或 水 接触 , 就 会 产生 
的 重要 安全 隐患 之 一 。KiKuChir4g 等 人 使 
能 高 于 低能 爸 反 应 路 径 
的 解 离 能 垒 受 到 水 分 子 相对 于 液体 钠 表 面 的 方向 的 影响 ， 
入 2at% 的 Ti 纳米 颗粒 来 抑制 
阐释 了 抑制 SWR 的 机 理 ， 即 由 于 钠 原 子 与 钛 纳米 颗粒 表 
水 合 能 量 和 启动 SWR 所 需 的 活 


tS REA BUR © Kim [等 人 发 现 可 以 在 液体 钠 加 


的 水 解 离 能 ， 


附 了 一 层 钠 原子 , 增加 了 


化 能 ， 且 减少 了 有 效 反 应 区 域 ， 从 而 抑制 了 钠 水 反应 。 


在 采用 金属 燃料 的 快 地 
是 一 种 设计 方案 ， 因 此 需要 考虑 金属 燃料 以 及 


钠 不 仅 与 氧气 、 水 蒸汽 发 生 
PBE 泛 函 ， 利 用 AMD 研究 了 液态 钥 


剧烈 反应 


， 核 燃料 
中 含有 镁 系 元 素 原子 的 三 


时 进行 了 模拟 计生 


分 别 为 -387.1 kJ-mol! 和 -374.0 kJ-mol!, 600 K 和 1 000 K 温度 下 Naz0 WRTA 28.6 


具有 较 好 的 一 致 性 ， 


剧烈 反应 。 钠 水 反应 (Sodium-Water Reaction, SWR) 
AIMD 对 SWR 的 初始 过 和 


SWR 


SECUS E 


裂变 元 素 也 会 对 其 性 能 产生 影响 。 


Li F5 


体 


二 元 金属 


性 质 ， 发 现在 723 K P. fic h3U pr SUB EXCESS. AER ANA 


中 的 扩散 率 对 温度 的 依赖 性 最 强 。Samin P? 


在 此 工作 基础 上 使 


SAS ABET SATE RT Fe. Th 


ot 
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| GGA- PBE iz if 
果 表 明 : Æ 1000 KIN, AEG A 


发 进行 。 计 算 
水 分 子 的 氢 原 


系统 中 ， 为 了 提高 燃料 元 件 内 的 热传导 ， 燃 料 与 包 壳 的 间隙 充填 液态 钠 


然 料 的 裂变 产物 与 钠 的 相 容 性 。 FERN, WAS 


1 GGA- 


AE 723~1 000 K FW 


, 利 


考 则 表现 为 液态 ， 钟 在 钠 
| AIMD 对 钠 - 
的 溶解 度 小 于 0.78 


at.%, FRRO 5.57x10? ms. DAA, b SCIT el RT E S SO RAS AT A 


18 hth Sophie HERRI DFT 计算 


TEENA REP, “AAG LBE Befthit zx RAR AB 
FI, A T SEU E E OS LA 
用 GGA-PBE 泛 函 进行 DFT EHI, MA T EURIBOR T TE Fe(111) 
别 为 一 
能 分 别 为 5.11 eV 和 4.58 (4.56) eV。Song50 等 人 使 
目 ， 铅 和 锐 原 子 在 (100) 表 面 


Js 


IU Ji T^ HS] 
表面 的 逸 上 


Bx 


7L 


T LBE 与 Fe(100) 表 面 的 相互 作 月 


能 分 


附 效果 最 好 ， 吸 附 


人 使 月 
TIE h BETA 
配 位 数 降低 所 


原子 在 Fe(111) 粗 糙 表 面 的 吸附 能 , 增加 


H GGA 泛 函 进行 DFT 模拟 ， 系 统 地 讨 


HB 


明 粗 糙 的 Fe(111) 表 面 比 至 
NM. mn $a. Uu 


it att, WETA 


的 液态 金属 腐蚀 。 除 了 进行 大 量 的 实 
央 ， 许 多 学 者 也 开展 了 模拟 计算 研究 工作 。Long59 等 人 使 

表面 的 相互 作用 ， 发 现 fec 位 点 对 铅 
2.66 eV 和 -3.07 eV。 铁 原子 在 洁净 表面 和 吸附 有 铅 〈 铅 ) 
| GGA-PBE 泛 函 进行 DFT 模拟 ， 研 
的 间 隐 位 置 都 有 良好 的 吸附 。Li523 等 
用 ， 吸 附 能 
面 的 铁 原子 


icf LBE 原子 与 
面具 有 更 强 的 ) 


Fe(111) 粗 糙 表 面 的 相互 作 
性 ， 这 主要 是 由 于 粗糙 表 
附 能 ， 发 现 增加 铝 含 量 可 显著 提高 铅 


AK 


f 


RUE Eu] Eque 


HÁ 
E Fe(111) 面 的 吸附 能 。 


Li 


iu53] 等 人 使 
RUBIA E IR a 


, 


] GGA 泛 函 进行 DFT 模拟 
界 渗透 导致 其 膨胀 ， 从 而 肖 


LiuP425 A fii 


腐蚀 的 敏感 性 。 


用 GGA 泛 函 进行 DFT 模拟 , 又 进一步 通过 吸 PB 
结果 表明 LBE 在 Fe 或 Fe304 表 面 


更 耐 LBE 溶 
位 置 被 稳定 吸 
锈 钢 
扩散 系数 比 前 


DingP925 A fii Hj GGA-PBE 泛 函 进行 DFT 模拟 , 研究 了 不 锈 钢 (110) 晶 面 


PE BS ah 


附 。Zhou6535 等 人 使 


J 合金 元 素 (Cr. Ni. Al. Si) 在 Fe(111) 


人 的 研究 要 低 得 多 ， 说 明 在 


1 弱 唱 界 内 聚 力 造 成 脆 化 。 不 锈 


。 此 外 ，LBE 原子 可 以 很 容易 地 吸附 在 Fe 表面 ， 但 
| GGA-PBE 泛 函 进行 DFT 模拟 ， 研 究 了 添加 O 原子 的 LBE 与 不 


粗糙 表 
研究 了 LBE 在 bee Fe25(013) 章 界 处 的 腐蚀 机 理 ， 发 现 铅 
表面 会 形成 Fe304 保护 层 ， 
和 逃逸 能 描述 了 Fe 和 Fe304 表面 对 LBE 
的 吸附 降低 了 Fe 的 逸 出 能 ，Fe304 表面 比 Fe zT 
能 在 Fe304 表面 上 的 某 些 


Hf 


L 
Li 


= 


而 的 界面 腐蚀 行为 ， 得 到 的 Pb、Bi 和 0O 原子 的 
面 处 形成 的 氧化 层 可 以 有 效 地 提高 钢 的 耐 蚀 性 。 
合金 元 素 和 氧 在 液态 


AX 
Fe-LBE 界 


LBE 中 的 扩散 行为 ,计算 出 Fe-O、Cr-O、Ni-O、Al-O 和 Si-O 键 的 形成 能 分 别 为 -0.56 eV. —0.60 
eV、0.05 eV.—1.08 eV 和 -1.11 eV, 并 发 现在 500 CCH AIMD 模拟 中 ,20ps 内 Ni-O 键 就 会 断裂 ，Fe-O， 
Cr-O, AI-O 和 Si-O 键 则 相当 稳定 。 最 终 解释 不 锈 钢 在 LBE 中 的 行为 。 钊 更 倾向 于 溶解 到 LBE 中 ， 
而 不 可 能 形成 镍 - 氧 络 合 物 。 铁 、 铬 、 铝 和 硅 可 吸附 在 钢 -LBE 界面 ， 并 可 与 氧 原子 结合 ， 形 成 氧化 
JZ. LBE 中 的 氧 原子 倾向 于 与 钢 -LBE 界面 区 存在 的 铝 和 硅 结合 。 因 此 ， 在 钢 中 挨 杂 铝 / 硅 有 助 于 形 
成 保护 性 氧化 层 ， 从 而 提高 钢 在 液态 LBE 中 的 耐 腐蚀 性 。 


19 ” 锂 冷 环境 中 材料 腐蚀 的 DFT 计算 


聚变 堆 第 一 壁 材料 主要 为 316 不 锈 钢 ， 流 动 液体 名 


良好 的 慢 化 


子 及 防止 


A RC RE ZH 
进行 


Sb 恒 
i H 


能 量 上 都 倾向 


| DFT BU, XI 


件 。 然 而 与 结构 材料 接触 时 
ET XE TCR IE REO IC D 


子 泄露 的 作用 ， 可 以 减少 燃料 


1 和 铁 原子 的 逃逸 进行 计生 


CERO Ay 


EE BLSE TS — BEAT BL AA AST BL, FUB 
循环 和 杂质 产生 ， 并 通过 提供 自 愈合 表面 来 
带 来 溶解 腐蚀 。Xu57 等 人 使 用 GGA-PBE iz ER 

。 结 果 表 明 ， 锂 原子 和 铅 原 子 在 


， 液 体 锂 会 


于 吸附 在 bec-Fe 表面 ， 并 进 


的 顺序 ， 表 了 
工作 环境 的 温 
bcc-Fe 中 的 扩 
<111>W 
0.063 eV. 


铃 位 点 的 形成 


铁 原子 逃逸 概 率 增加 。 铁 的 
度 有 很 强 的 依赖 性 。Li5 使 


Ve AES OT eT A. BR A A 


步 加 速 了 表面 铁 原 子 的 逃逸 。 按 (110)<(110)<(001)<(111) 
或 铅 原子 的 覆盖 范围 和 
| GGA-PBE 泛 函 进行 DFT 模拟 , 进一步 研究 了 锂 原子 在 


散 行为 ， 发 现 取代 位 是 锂 原 


S 


子 的 首选 位 置 ， 通 过 空位 机 制 扩 散 。 在 所 有 的 间隙 位 置 


, 


低 ? 间 际 锂 Jf 


Sb E, 


He HX 
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f 


的 扩散 是 通过 其 在 <111> 位 点 的 旋转 来 完成 的 ， 能 垒 仅 有 


反应 堆 合金 材料 在 熔 盐 堆 环境 下 腐蚀 的 DFT 模拟 计算 


第 四 代 反 应 堆 中 ， 熔 盐 堆 在 针 资 源 利 用 和 核 废料 处 理 方面 具有 显著 的 优势 ， 并 因 其 固有 安全 性 
高 、 防 核 扩散 等 优点 ， 从 而 受到 广泛 关注 。 其 所 使 用 的 炊 盐 主要 有 氟 化 盐 与 氧化 盐 。 以 Hastelloy N 
合金 为 代表 的 镍 基 合 金 ， 由 于 其 良好 的 高 温 机 械 性 能 和 优异 的 耐 腐蚀 性 ， 在 炊 盐 堆 中 被 用 作 主 要 结 
构 材 料 。 熔 盐 堆 环境 下 反应 堆 合 金 材料 腐蚀 主要 是 裂变 产物 〈 主 要 是 确 ) AEC RMT MURA 
金 元 素 向 熔 盐 的 扩散 。 表 4 汇总 了 熔 盐 堆 环境 下 反应 堆 合 金 材料 腐蚀 的 DFT 文献 情况 。 
表 4 熔 盐 环境 下 反应 堆 合 金 材料 腐蚀 DFT 计算 文献 汇总 


Table 4 Literature on DFT calculations of corrosion of reactor alloy material in molten-salt environment 


研究 机 构 第 一 作者 主要 内 容 年 份 
Research institution First author Research content Year 


Shanghai Institute of Applied W.G.Liu T t RED & X SODAR RARI) 2014 
Physics JE RIT a REDI 
The process of Stress corrosion cracking of 


tellurium-containing Ni-base alloy © 5(012) 
symmetric inclined grain boundary 


Shanghai Institute of Applied W.G.Liu HT SHEN Ni(111) 表 面 和 5(012) 晶 界 的 应 力 2014 
Physics 腐蚀 开裂 过 程 的 影响 [0 
Effect of tellurium on the stress corrosion 


cracking process of niobium-containing Ni(111) 
surface and 225(012) grain boundary 


Shanghai Institute of Applied W.G.Liu A RA TG Be Bite IML 73 JS ORT AREA] RR] IU 2014 
Physics Effects of various elements on stress corrosion 


cracking induced by tellurium 


Shanghai Institute of Applied L.Mi Titi PEAR ZEB SE i FEAL AY OD A EF A 2019 
Physics Distribution and existence forms of tellurium at 
the grain boundaries of nickel-based alloys 


Shanghai Institute of Applied H. Ai Ni-26W-6Cr 合金 和 GH3535 合金 在 800°C FYE — 2019 
Physics FLiNaK 熔 盐 中 的 腐蚀 行为 [9] Corrosion 


behavior of Ni-26W-6Cr alloy and GH3535 alloy 
in FLiNaK molten salt at 800°C 


University of California N. Winner 三 种 氟 化 物 在 添加 和 不 添加 Cr 时 的 结构 的 2021 
Structure of three fluorides with and without Cr 
Massachusetts Institute of QJ.Li NaCl - CrCl; If] £& 44-5; 4 ^f: HER! 2021 
Technology Structure and chemical properties of NaCl—CrCl 
Georgia Institute of K.Hanson RS KCI-MgCbL JE di zh AY fc vl Reaction 2022 
Technology of nickel with KCl-MgCl eutectic salt 
Wright-Patterson Air Force D. Tyler NizoW»oNbio TE SCELTA HE EI] 5 89] VR RT 70071 2023 
Base Early dissolution behavior of NizoW»oNbio in 


chloride molten salt 
EE BAS Ae BIAS ERE ER AN, Hastelloy N 合金 可 以 观察 到 明显 的 品 间 裂纹 。 这 些 裂纹 的 


Ui BUE ET E ih I MIRESE 


体 


Li 


= 


u 


o LiuP? 
Tii FETE LYSOL) MRA FE. OR FE Mi S| ERE 
先 占据 替代 位 ， 由 于 原子 尺寸 的 不 匹配 ， 磋 引起 了 唱 界 的 膨胀 ， 从 
0 等 人 使 用 GGA-PBE 泛 函 进行 DFT 模拟 ， 研 究 了 锟 掺 杂 后 的 Ni(111) 表 面 和 也 
包 倾 向 于 向 晶 界 偏 析 并 形成 较 强 的 Nb-Ni 键 , 阻止 Te 向 合金 内 部 迁移 从 而 
金 对 确 致 晶 间 脆 化 的 抵抗 能 力 。 在 此 基础 上 ，Liut61 使 月 


等 人 使 


的 应 力 腐蚀 


开裂 。 


现 稀土 元 素 与 磷 元 素 同 样 能 


mr Z8 BE 5(012) dà 


在 铝 基 合金 唱 界 处 的 空间 分 布 和 化 学 形态 。 


FT BEXT 


而 削弱 了 界面 的 Ni-Ni 键 。 


| GGA-PBE 泛 函 进行 DFT 模拟 , 通过 模拟 
证 实 了 确 在 唱 界 区 域 的 优 


随后 ， 


上 75S(012) 唱 界 ， 发 现 


u-XRF E 


后 的 u-XRD 分 析 证 实 了 唱 界 腐蚀 
电荷 从 磋 的 2p 轨道 转移 到 镍 的 4s 轨道 ， 
反应 过 程 如 下 所 示 [9， 


产物 为 确 化 镍 (NisTez) 和 镍 - 
DFT it & 


hii 


也 表 


BH tit El 


l 705°C |. 
6NiTe (c) ——> 2Ni3Tez (c) + Te» (g) 


. 800€ 
NisTe2 (c) ——93Ni (c) + 2Te (g) 


+ 128 


由 于 活性 合金 元 素 的 选择 性 溶解 ， 结 构 合 金 暴 露 在 熔 盐 
有 化 学 和 使 用 GGA-PBE 泛 函 进行 DFT 模拟 ， 研 究 了 Ni-26W-6Cr 合金 和 GH3535 合金 


方法 ， 结 合 : 


.9 kJ-mol! 


在 800 °C FE FLiNaK 盐 中 的 腐蚀 行为 。 


结果 表明 ， 


增强 上 
H GGA-PBE 泛 函 进行 DFT 模拟 ， 进 一 步 发 
唱 界 的 脆 化 /强化 作用 与 偏 析 谈 度 密切 相关 ， 
而 铬 没有 明显 的 偏 析 趋 势 。 随 后 Mi 等 人 使 用 GGA-PBE iz PRETI DFT 模拟 ， 结 合 实验 研究 了 确 
谱 显 示 ， 确 的 浓度 从 唱 界 向 晶 粒 逐渐 降低 。 


会 受到 腐蚀 。AilGl 等 人 采用 项 浸 腐 蚀 


nat, dern RAE 


随 


回 溶 体 。 特 别 是 ，u-XANES 结果 显示 
YES rp BS ga-t XR Ae sg 


价 键 。 


两 种 合金 的 表面 均 存 在 Cr 损耗 层 ， 最 大 深度 约 


为 13 um, W 和 Mo 的 损耗 非常 轻微 .Winner 等 人 [的 使 用 GGA-PBE 泛 函 ,利用 AIMD 研究 了 2KF-NaF、 
2LiF-BeF; 和 3LiF-AlF; 三 种 气 化 物 在 添加 和 不 添加 Cr 时 的 结构 ,他们 发 现 Cr 可 以 被 氟 离 子 溶剂 化 ， 
Cro、Cr、Crn 在 三 种 氟 化 物 中 都 可 以 被 不 同 数量 的 F 配 位 。2KF-NaF 熔 体 表现 出 高 度 的 短程 有 序 ， 


而 2LiF-BeF» 和 3LiF-AIF; 中 程 有 序 。,Lits5 等 人 使 月 


的 结构 与 化 学 性 质 ， 发 ] 
从 而 直接 影响 腐蚀 速率 。 
除了 铬 元 素 之 外 ， 


其 他 合金 元 素 也 会 影 


究 了 钊 与 KCI + MgCh 共 唱 盐 混合 物 的 反应 ， 发 现 Mg? 57 


响 镍 基 合 金 与 熔 盐 的 反应 。Hansonl66 
C PA MgCE ITE SN ES ERA TR 


的 熔 盐 环境 
BERT ULL S SEO E 1 rh 
在 氧化 盐 下 的 早期 溶解 行为 。 发 现 
从 而 增强 了 镍 基 合 金 抵抗 熔 盐 府 包 


的 能 


结语 


DFT 已 经 在 反应 堆 合 金 材料 腐蚀 


， 握 原子 需要 与 其 他 熔 盐 离子 竞争 ， 
o Tylerl$7125: A fiti H 


岗 表 面 附近 的 Cr 浓度 最 高 ， 会 影响 盐 的 扩散 系数 、 粘 度 和 氧化 还 原 


H GGA-PBE 泛 函 进行 DFT 模拟 ,研究 了 NaCl-CrCl 


HA, 


等 人 使 用 AIMD 研 


X 


He AS XE 


领域 得 到 了 一 些 应 月 


ESTNE 
[8] fs HH 


LY 


容易 与 氧 结合 ， 并 且 


TET 


结合 。 在 真正 


因此 ， 向 熔 盐 中 添加 过 量 的 
H GGA 泛 函 ,研究 了 钊 高 温 合金 NizoW2oNblo(110) 


形成 的 化 合 物 不 容易 浴 解 于 氧化 盐 ， 


,主要 包括 : 水 冷 堆 环 境 下 铬 合金 的 氧化 、 
致 应 力 腐 蚀 与 其 他 材料 的 腐蚀 ; 液态 金属 冷却 堆 环境 下 钠 -水 反应 、 铬 铅 溶 解 腐蚀 与 锂 原 


f 


合金 材料 的 扩散 ; 
DFT Æ RMH 


入 合金 材料 腐蚀 研究 中 主要 用 于 


模拟 : 腐蚀 


离 ， 腐 蚀 介质 元 素 在 反应 堆 合金 材料 表面 及 


Sb Ei, Lr. 
能 影响 等 o 


Sei, HATHI 
DFT 模拟 研究 依然 偏 少 ， 
合金 的 SCC、IASCC 机 理 ， 
下 腐蚀 性 能 预测 等 ， 其 


YR, DFT 模拟 通常 只 能 


题 大 多 都 发 生 在 高 温 条 件 下 ，DFT 模拟 难 


以 对 此 进行 精确 


原子 尺度 研究 ， 但 由 于 其 计生 
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晶体 内 部 的 扩散 行为 ， 腐 蚀 


环境 下 合金 元 素 向 炊 盐 的 扩散 与 裂变 产物 确 引 起 的 应 力 腐蚀 


介质 元 素 在 反应 堆 合金 材料 表面 的 吸 


开裂。 
附 与 解 


合金 材料 腐蚀 的 DFT 模拟 计算 还 存在 一 些 问 题 。 首 


些 腐蚀 问题 仍然 需要 开展 DFT 模拟 计算 ， 如 水 冷 夫 
高 温 气 冷 堆 环境 下 结构 材料 的 氧化 机 理 ， 高 箭 合 


介质 对 反应 其 


H 


合金 材料 的 性 


J 


该 领域 已 开展 的 


"m 


EASE RAG. PRAE 


金 


等 新 材料 在 堆 内 环境 


研究 OK 温度 下 的 情况 ， 反 应 夫 


环境 所 涉及 的 腐蚀 问 


kh 


时 间 较 长 且 成 本 相对 较 高 ， 往 往 被 局 FP 


究 ; Hoa,» DFT 模拟 虽然 可 以 开展 电子 、 
R 在 较 小 的 时 间 、 空 间 尺 度 中 ， 


O 


N 
PA 


能 选择 简化 后 的 典型 模型 来 开展 特定 过 程 的 腐蚀 计算 ， 模 拟 过 程 中 考虑 的 模型 与 反应 介质 较为 单 


* 
反应 堆 材料 腐蚀 问题 将 具有 更 为 广阔 的 空间 。 未 来 还 有 以 下 几 个 方面 工作 待 开展 : D 进一步 扩大 反 


mi 


难以 充分 揭示 复杂 的 材料 腐蚀 问题 ， 所 得 结果 也 不 一 定 能 够 与 宏观 腐蚀 试验 结果 相符 合 。 
随 着 计算 机 技术 和 计算 方法 的 发 展 ，DFT 计算 的 精度 与 速度 都 在 不 断 提 高 ， 利 用 DFT 计算 和 


人 材料 腐蚀 问题 DFT 模拟 的 研究 范围 ， 如 上 述 堆 型 环境 其 他 腐蚀 问题 、 其 他 堆 型 环境 腐蚀 问题 、 
新 材料 选材 论证 等 ， 2) 基于 材料 基因 工程 中 理论 计算 、 实 验 、 数 据 库 协同 发 展 的 理念 ， 结 合 机 器 学 


习 与 高 通 量 计算 方法 ， 开 展 反 应 堆 材 料 腐蚀 问题 研究 ，3) 针对 DFT 计算 尺度 小 的 问题 ， 结 合 其 他 
模拟 方法 (分 子 动力 学 、 相 场 模 拟 和 有 限 元 方法 等 ) 开展 反应 堆 材 料 腐蚀 的 多 尺度 模拟 计算 。 
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